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FOREWORD 

Continuing  Interest  In  the  possibility  of  using  clustered ^varhesd 
weapons  In  non-nuclear  warfare  has  given  sustained  popularity  to 
NAVORD  REPORT  7019  (A  Handbook  on  the  Effectiveness  of  Cluster  Weapons 
Against  Unitary  Targets,  by  Eldon  L.  Dunn).  Convenient  methods  of 
extending  the  results  to  larger  nuBft>ers  of  bcslblets  and  larger  dispersions 
without  additional  machine  casiputatlon  have  been  employed  recently  at  the 
U.  S.  Naval  Ordnance  Test  Station.  It  Is  believed  that  these  methods  will 
be  of  Interest  to  \isers  of  the  "Handbook”  and  others  faced  with  the  prob¬ 
lem  of  estimating  the  effectiveness  of  cluster  weapons. 
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INTRODUCTION 

Interest  in  the  possibility  of  increasing  the  effectiveness  of  con¬ 
ventional  weapons  by  employing  clusters  or  salvos  of  many  borablets  has 
been  accompanied  by  corresponding  interest  in  the  problem  of  estimating 
the  over-all  kill  probability  of  cluster  weapons.  An  effective  approach 
to  this  problem  has  appeared  in  NAVORD  REPORT  7019  by  Eldon  L.  Dunn,l  in 
which  the  results  of  a  parametric  study  of  cluster-weapon  kill  prob¬ 
ability  against  unitary  targets  are  presented  in  handbook  form,  A  uni¬ 
tary  target  is  one  that  can  be  said  to  be  destroyed  (or  damaged  to  a 
specified  extent)  as  a  unit;  it  is  to  be  distinguished  from  a  target  or 
group  of  targets  that  can  be  said  to  be  damaged  over  a  fraction  of  its 
total  area  (area  target).  Only  unitary  targets  are  considered  iri  this 
report. 

Although  the  graphs  in  the  Handbook  (NAVORD  REPORT  7019)  are  quite 
useful  in  estimating  the  effectiveness  of  proposed  cluster  weapons,  their 
'applicability  is  sometimes  limited  by  the  difficulties  of  interpolation 
and  by  the  absence  of  data  for  more  than  kOO  bomblets  or  for  delivery 
standard  deviations  greater  than  8  target  widths.  It  has  been  observed, 
however,  that  some  simple  relationships  appear  to  exist  among  the  para¬ 
meters  in  the  Handbook  within  certain  ranges  of  their  values.  Simple 
rules  can  be  derived  for  extending  the  ranges  of  some  variables  by  using 
"artificial"  values  of  other  variables  with  the  graphs  provided  in  the 
Handbook. 

In  this  report,  these  empirically  determined  rules  are  briefly  de¬ 
scribed  and  Illustrated.  Graphs  based  on  these  rules  are  presented  to 
permit  salvo  kill  probabilities  to  be  estimated  without  using  the  Hand¬ 
book,  for  many  cases  both  within  and  beyond  the  scope  of  the  Handbook 
graphs . 

Ihen,  in  an  attempt  to  provide  insight  into  the  basic  relationships 
and  to  verify  their  reliability  in  extrapolating  the  Handbook  data,  the 
basic  equation  used  in  the  preparation  of  the  Handbook  is  examined.  Two 
simple  approximations  and  some  changes  in  notation  are  Introduced  to  aid 


^  U,  S,  Naval  Ordnance  Test  Station.  A  Handbook  on  the  Effective¬ 
ness  of  Cluster  Weapons  Against  Iftiitary  Ihrgets,  by  Eldon  L,  Dunn, 
China  Lake,  California,  NOTS,  12  January  I960.  (NAVORD  REPORT  7019, 
NOTS  TP  2382^ 
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the  discussion,  which  generally  supports  the  application  of  the  empiri¬ 
cal  rules  to  cases  where  standard  deviations  and  numbers  of  bomblets 
are  relatively  large.  It  appears  that  the  basic  equation  can,  for  such 
cases,  be  "condensed"  into  a  form  in  which  the  salvo  kill  probability 
is  approximately  a  function  of  Just  two  variables  that  correspond  to 
those  that  appear  in  the  empirically  derived  graphs. 

In  an  appendix,  a  few  minor  errors  in  the  Handbook  are  briefly 
listed  for  its  users*  convenience. 


THE  EMPIRICALLY  DETERMINED  RULES 

The  observations  described  in  this  section  essentially  make  up  the 
foundation  of  this  report.  While  the  mathematical  discussion  is  Impor¬ 
tant  in  understanding  and  generalizing  the  observed  relationships,  both 
the  original  inspiration  and  the  detenoination  of  approximate  limits  of 
application  have  come  from  Inspection  of  the  Hemdbook  graphs.  The 
reader  who  has  a  copy  of  the  Handbook  and  is  Interested  only  in  some 
simple  methods  of  Interpolating  and  extrapolating  the  data  may  find  this 
section  adequate  for  his  needs.  On  the  other  hauid,  the  reader  vdio  does 
not  have  at  heuid  a  copy  of  the  Handbook  is  encouraged  to  pass  quickly 
over  the  many  references  to  the  Handbook,  since  the  graphs  and  equations 
that  appear  later  in  this  report  can  be  used  independently. 

The  problem  treated  in  the  Handbook  is  that  of  finding  the  salvo 
kill  probability  of  a  cluster  weapon  whose  bomblets  are  normally  dis¬ 
tributed  about  a  point  that  is  itself  normally  distributed  about  the 
aimpoint,  against  a  rectangular,  unitary  target. 

The  graphs  in  Appendix  A  of  the  Handbook  show  the  salvo  kill  prob¬ 
ability  ( /*  )  as  a  function  of  the  ammunition  standard  deviation  (  ) 

for  various  values  of  delivery  standard  deviation  (  ),  values  of 

probability  that  a  hit  by  a  bomblet  is  a  kill  i  /C/r),  numbers  of  bomblets 
(  A/  ),  and  target  lengths  (  L  ).  The  values  of  ,  and  L 

are  given  in  terms  of  target  width  (  ).  (Specifying  the  x  -direction 

standard  deviations,  and  $  ^^^o  specifies  the  /  -direction 
quantities,  since  the  ratios  fixed  at  1.^ 

in  the  Handbook. )  The  values  of  that  correspond  to  maximum  values 
of  P  on  each  curve  are  referred  to  here  as  optimua  values.  The  graphs 
in  Appendices  C  and  D  of  the  Handbook  show  Just  the  maximvoi  salvo  kill 
probability  )  as  a  function  of  At  or  of  ^  ,  when  is  always 
optimum. 

Certain  similarities  among  curves  frm  these  Handbook  graphs  become 
apparent  when  they  are  used  extensively.^  For  exanq>le,  if  V  is  changed 


^  These  similarities  were  first  pointed  out  to  the  author  by  M.  L. 
Bralthwaite  of  the  Weapons  Development  Department  at  the  U.  S.  Naval 
Ordnance  Test  Station. 
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but  Kfi  is  also  changed  so  that  the  product  Is  not  changed^  then 
the  value  of  P  often  is  not  changed.  Thus,  the  curve  for  N  ■  200  and 
Kfi  m  0.4  from  Fig.  A59  of  the  Handbook  is  about  the  same  as  the  curve 
for  A/  ■  400  and  Kp  •  0.2  from  Fig.  A60.  Further,  the  same  curve  will 
hold  for  A/  ■  800  and  Kp  ■  0.1  or  for  AJ  m  0.05,  if  target 

length  remain  2^  and  3 IV, respectively.  (This  has  been  verified  by 
separate  machine  computation  with  a  000  and  l600. )  It  appears,  then, 
that  the  sin^lest  method  for  extrapolating  the  Handbook  data  to  larger 
values  of  N  is  to  replace  the  desired  values  of  N  and  Kp  with 
"artificial"  values  hj*  and  Kp*  that  appear  in  the  Handbook  and  satisfy 

N*Kp*  -  NKp . 

Example  1:  Find  the  maximum  probability  {  P^)  for  a  target 
20  X  60  ft  with  ffpx  »  80  ft,  -  120  ft, 

A>  «  0.025,  and  A/-  l600,  3200,  ^00. 

Since  1  /c^  ■  1.5»  and  the  target 

length  is  IW  $  Fig.  Cll  in  the  Handbook  is  used  with 
Kp*  «  0.4  and  N*»  100,  200,  400.  Then  »  0.35,  0.52, 
0.70. 

Relationships  involving  other  parameters  are  also  observed,  but  they 
do  not  hold  as  generally  as  the  one  that  involves  just  /V  and  Kp,  If 
large,  it  can  be  seen  that  P  is  approximately  the  same  if  the 
product  A/Kp  L  is  the  same,  where  L  is  the  target  length.  Compare, 
for  example,  the  curve  for  ■  0.6  in  Fig.  A29  of  the  Handbook  with 
the  curve  for  ■  0.1  in  Fig.  A66.  In  each,  case,  m  k  IV  and  the 

product  A/Kp  L  »  120.  The  curves  are  not  identical,  but  the  values  of 

p  differ  by  only  about  0.01.  Under  conditions  when  is  large 

enough,  then,  artificial  values  A/*,  K^ ,  L*  that  can  be  found  in  the 
Handbook  graphs  can  be  substituted  for  desired  values  ^  Kpj  L  it 
A/*  Kp^L*  - 

Similarities  among  curves  with  different  values  of  also  be 

observed.  For  example.  Fig,  A28  for  /V  ■  100  appears  to 

be  nearly  the  same  as  Fig,  A36  for  •  8lVand/V»  400,  It  is  necessary, 

however,  to  change  the  scales  to  scales  of  ,  80  that  all 

graphs  in  Appendix  A  of  the  Handbook  have  the  same  horizontal  scale. 

Since  a  cheuige  by  some  factor  in  often  is  equivalent  to  a  change 

by  the  square  of  the  factor  in  AA  or  ,  it  appears  that  in  many  cases 
P  is  approximately  the  same  if  the  quantities  NKpL/(^^  and 
are  the  same.  The  curves  for  which  NKpL/p^^*’  is  the  same  are 
nearly  eU.ike  if  is  large  enough,  but  for  smaller  the  curves 

begin  to  diverge  for  the  smaller  values  of  Cjjgy  /^x  •  ^  general, 
the  values  of  ^  agree  within  about  0,02  if  large  enough  that 

either  is  not  less  than  Its  optimum  value  or 

larger  than  about  half  the  corresponding  target  dimensions.  Also,  A/ 
should  be  larger  than  about  10, 
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Example  2:  Find  P  for  a  target  20  x  60  ft  with  ■  160  ft, 

-  0.6,  and  V  ■  l600,  if  crjk^  -  ^0  ft,  240  ft 
(  m  1.5),  Also  find  the  opti¬ 

mum  ^lue^of  ^^■^auid  the  corresponding  maximum  prob¬ 
ability,  ^  ,  It  is  seen  that  ■  8^.  By  adopting 

the  artificial  values  N* m  400  and 
Fig.  A66  in  the  handbook  can  be  used.  ]rar  a 

(that  is.  4o  ft  or  IM/),  the  curve  i^r  Kp  a 
0.6  gives  P  a  0.34;  for  a  1.5  (that  ia,  240  ft 

or  6t^XP"0,7^»  The  optimum  value  of  is  about 

1.05  or  about  l68  ft  and  the  maximum  probability 

is  about  0.79*  From  separate  machine  computation,  the 
results  are  Pm  0.32  and  0.74  for  m  40  ft  and 

240  ft,  and  a  0.795  at  »  165  ft. 

If  only  the  maximum  values,  ^  ,  are  desired,  the  graphs  in  Appendices 
C  and  D  of  the  Handbook  are  used  with  any  combination  of  artificial 
values  that  satisfy  ^  -  NKpL/O^^^  .  (The  ratio 

i'ts  optimum  value  in  this  case.)  Ihe  next  two 
exaiiQ)les  snow  how  easily  some  of  the  examples  in  the  Hwdbook  are  worked 
by  this  method.  In  both  cases,  *  1.5  • 

Example  3!  Find  for  a  teurget  4o  x  120  ft  with  a  60  ft 

(a  1.5/t/),  Kp»  0.5#  and  N  »  10.  Since  curves  for 
Kp  m  0.5  do  not  appear  in  the  Handbook,  Kp  a  0.6  is 
used  with  a  1.641V.  The  value  of  P^  is  then  0.35# 

read  from  Fig.  06  of  the  Handbook.  (The  linear  inter¬ 
polation  suggested  on  page  7  of  the  Handbook  gives  0.34 
for  . ) 

Example  4;  Find  for  a  target  10  x  20  ft  with  a  30  ft 

(■  3  ^)#  Kp  a  0.45,  euid  N  a  120.  In  this  caise,  none 
of  the  parameters  has  a  value  that  appears  in  the  Hand¬ 
book.  Letting  Kjt  ■  0.4.  L*  a  3,  and  ^  a  kw, 
is  coBiputed  from  ^*(0.4)(3)/l6  a  (l20)(0.45)(2)/9. 

Then  N*  a  16O,  and  a  0.46  from  Fig.  Oil  of  the  Hand* 
book.  A  separate  machine  conqputation  also  gives  0.46 
for  /2  . 

It  is  easily  seen  that  this  method  is  less  conqplicated  and  more  ac¬ 
curate  than  the  interpolation  procedure  by  which  the  same  problem  as 
Example  4  is  attacked  on  pp.  7-9  of  the  Handbook. 

If  the  ratio  ^y/<7p^  is  not  1.5#  the  formula  for  artificial  values 
can  be  written  in  the  form 

(1)  = 
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vnien  finding  values  of  ^  extension  of  the  "sin^lifled 

method"  described  on  pp.  9*13  of  the  Handbook.  If  A/*  ■  A!  and  ■ 

Kp  t  (l)  above  is  the  same  as  Eq,.  (2)  on  p.  11  of  the  Handbook,  since  K 
in  the  Handbook  is  and  esc  S  in  the  Handbook  is  •  Of 

course,  vdien  finding  P  from  the  curves  in  Appendix  A  of  the  Handbook, 
the  condition 

(2 )  observed, 

Exanq>le  5:  Find  P«  for  the  same  conditions  as  Exaaqple  4,  except 

A/** is  2^0  and  the  curve  for 

Kp»  0.4  from  Pig,  Cll  of  the  Handbook  gives  ■  0.57» 
This  is  sklso  the  correct  (machine-con^uted)  value. 
Similarly,  if  ■  2,  N*  is  the  same  as  the 

original  A/,  120,^  and  from  the  same  curve  is  0.39. 

Values  of  P^  from  the  Handbook  graphs  are  accurate  within  about  0.01 
if  N  is  greater  than  10,  is  greater  than  ^tv',  and  dp  is  greater 

than  Values  of  P  are  generally  accurate  within  about  0.02  if  N 

is  greater  than  10  and  either  ^®®®  than  optimum  or 

d^^  and  dty  are  greater  than  i&out  naif  the  corresponding  target 
dimensions. 


METHODS  FOR  WHICH  THE  HANDBOOK  IS  NOT  REQUIRED 

Ibe  methods  Just  discussed  for  interpolation  and  extrapolation  of  the 
Handbook  graphs  suggested  the  possibility  of  "condensing"  some  of  the 
data  of  the  Handbook  into  a  few  curves.  Pursuing  this  idea,  if  is 
approximately  the  same  when  the  quantity  A/ATp  L / is  ^  same, 
it  seems  natural  to  adopt  a  single  new  parameter  Z  ,  where 

(3) 

and  plot  a  single  curve  showing  Pp  ae  a  function  of  Z.  Such  a  curve 
appears  here  in  Fig.  1.  Values  from  the  Handbook  curves  lie  very  close 
to  this  curve  if  A/  is  larger  than  about  10  and  is  larger  than  both 

M^/2  and  L/3,  In  general,  of  course,  since  the  iSilo  1»5 

is  included  in  the  Handbook  data,  the  equation  for  Z  should  Include  both 
d^  and  dp^  .  Also,  the  vailues  of  ,  dp^  ,  and  L  need  not  be 
expressed  in  target  widths  (iV)  if  lA/  appears  in  the  numerator  of  the 
foraula.  The  product  LW  is  Just  the  target  area,  T ,  so  the  general 
form  of  (3)  is 

(4)  Z  -  /.S  NKp  T /d^ 

When  ^  is  read  from  Fig.  1  using  Z  computed  from  (4),  it  is  accurate 
within  about  0.01  if  A'  is  greater  than  10  and  dp  and  dpy  are 
greater  than  about  half  the  corresponding  target  cumensions. 
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Examples:  For  tbe  previous  examples^  is  easily  found  from  Fig.  1 
as  follows: 

*  Ip  example  1.  for  A/  ■  l600,  3200,  and  6400,  the  values 
of  ^  are  7,5,  15,  and  30.  From  Fig,  1,  -*3  ■  0.35,  0.52, 
and  0.70,  in  agreement  with  the  values  from  the  Huidbook. 

.  In  example  2.  Z  ■  45,  and  from  Fig,  1  is  0,80,  ISie 
separate  machine  confutation  gave  0.793* 

.  In  example  3.  the  value  of  A/  is  rather  low  for  this 
metiiod;  still,  the  value  Z  ■6,67  gives  ■0.33  from 
Fig.  1.  value  from  Fig.  D6  of  the  Handbook  is  O.35. 

.  In  example  4,2«  12,  and  ^  ■  0,47  from  Fig,  1.  A 
separate  machine  confutation  gives  0.46. 

.  In  exaafle  5,  for  ■  1  and  2,  Z  ■  I8  and  9. 

The  values  of  ^  from  Fig,  1  are  0,57  and  0,39,  in 
agreement  with  the  previous  results. 

Although  Fig.  1  is  very  useful  for  estimating  values  of  ,  it  is 
not  applicable  if  other  than  optimum  vsilues  of  auid  cff  ,  the 

ammunition  standard  deviations,  are  ccmsidered.  Even  when  binding  /?  , 
it  does  not  indicate  vbaX  the  optimum  values  of  and  are.*  It 

appears,  however,  that  the  salvo  kill  probability,  P  ,  may  be  regarded 
in  many  cases  as  a  function  of  Just  two  variables,  Z  and  Af,  where  z 
is  defined  in  (4)  and  M  is  the  ratio  of  ammunition  standard  deviation 
to  delivery  standard  deviation.  That  is,  if 

P  can  be  plotted  as  a  function  of  A/f  for  various  values  of  Z •  Such  a 
fatally  of  curves  appears  here  as  Fig.  2,  In  general,  each  curve  is 
adapted  from  one  in  the  Handbook  for  which  Z  has  tbe  desired  value  and 
the  target  dimensions  are  small  with  respect  to  the  standard  deviations. 
Vedues  of  P  estimated  from  Fig.  2  are  accurate  within  about  O.03  if  N 
is  greater  than  10  and  either  M  is  not  less  than  optimum  or  and 

C\y  sure  greater  than  about  half  the  corresponding  target  dimension. 

The  curves  are  also  useful  for  estimating  optlmim  values  of  a/.  It 
appears  that  the  ammunition  standard  deviations  are  required  to  be  about 
the  same  size  aa  the  delivery  standard  deviations  (a/  is  near  l)  if  h^jg 
values  of  salvo  kill  probability  are  to  be  efficiently  achieved. 

Examples:  In  several  of  the  examples  treated  eeurller,  F;^  2  is  now 
used  to  estimate  optimm  values  of  emd  ,  aa 

follows: 
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.  In  exaaple  li  optimum  values  of  K1  for  z  ■  7.3>  15# 
and  30  are  approximately  0.63#  O.8O,  and  1.0^  respectively. 
Since  is  80  ft  >  the  optimven  values  of  are 

32  ft|  64  ft,  and  80  ft;  the  corresponding  optimiai  values 
of  are  78  ft,  96  ft,  and  120  ft. 

.  In  example  3»  where  Z  ■  6.67#  the  optimum  value  of  M 
is  about  0.62.  For  C7^  of  60  ft,  optimum  values  of 
and  d5e^  are  37  rt  and  36  ft,  respectively. 

.  In  examples  4  and  5.  the  optimum  ammunition  dispersions 
axe  easily  determined  in  the  same  manner; 


(Target  10 

X  20  ft. 

<r^ 

-  30  ft, 

Z 

fs 

a7" 

1.5 

12 

0.47 

1.0 

18 

0.37 

0.64 

2.0 

9 

0.39 

0.7 

A/  m  120,  /[;,  m  0.43) 


Optimum  V6Q,ues  of 


Ft 

ft 

23 

35 

23 

25 

21 

42 

It  has  been  assumed  here  that  the  ratio  of  to  is  edvays 

the  same  as  the  ratio  of  to  .  It  seems  reasonable  to  siq^se 

that  the  same  characteristics  of  an  attack  that  produce  an  elliptical 
distribution  of  delivery  error  will  produce  a  similarly  elliptical  dis¬ 
tribution  of  bomblet  inlets.  The  same  assumption  is  made  implicitly  in 
the  Handbook,  and  will  be  Introduced  again  in  the  next  section  to 
8iJiQ>llfy  the  mathematical  derivation.  It  is  this  ccunon  ratio,  M  ,  that 
is  optimized  when  the  maximum  probability  for  a  given  value  of  Z  is 
found.  This  ratio  (but  not  necessarily  the  separate  components  and 

)  is  also  optimized  in  the  Handbook  when  is  optimized. 


vrhile  the  formulas  and  figures  Just  presented  make  possible  the  esti¬ 
mation  of  values  of  salvo  kill  probability  without  further  use  of  high¬ 
speed  digital  computers,  and  ^ile  the  estimated  values  agree  well  with 
the  results  of  special  machine  computations  for  larger  numbers  of  bomb- 
lets  and  larger  standai^  deviations,  a  degree  of  mathematical  Justifi¬ 
cation  seems  desirable.  In  the  next  section,  the  bMic  equation  used 
in  the  preparation  of  the  Handbook  is  briefly  examined  in  the  lifi^t  of 
the  empirical  rules  and  methods  already  described. 


DISCUBSIOlf  OF  THE  BASIC  EQUATION 

As  described  on  pp.  2-4  of  the  Handbook,  the  problem  is  to  find  the 
probability,  P  ,  that  at  least  one  round  of  a  salvo  (cluster)  will  kill 
the  target.  If  the  aim  point  of  the  fire-control  system  and  the  center 
of  the  target  are  both  at  the  origin,  the  basic  equation  on  p.  2  of  the 
Handbook  may  be  written  as 
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Here,  PL  Is  the  probability  that  at  least  one  round  of  a  salvo  whose 
projected  impact  point  is  {  x,  y  )  will  kill  the  target.  That  is, 

(7)  PL  =  -  P  (Xj/)] 

where  p  {x  ,  y  )  is  the  probability  that  a  particular  round  of  the  sadvo 
will  kill  the  target.  Since  the  rounds  are  normally  distributed  about 
the  projected  impact  point, 

exp  (-iv^)cLy 

where  j 

g  ,  (.-x+-^XL)/oi^ 

c.(-y-^  rL)/o^^ 

As  in  the  Handbook,  XL  and  YL  are  the  dimensions  of  the  target,  //  is 
the  number  of  rounds  (bomblets),  Kp  Is  the  probability  that  a  hit  is  a 
kill,  (7p^  and  are  the  delivery  standard  deviations,  and 

euid  Opy  sre  the  amunition  (round)  stcuidard  deviations. 

It  is  con^ratlvely  easy  to  show  that  a  change  in  the  number  of 
rounds  can  be  approximately  ccmpensated  in  many  cases  by  an  inverse 
change  in  the  probability  that  a  hit  by  a  bomblet  is  a  kill.  If 
p{xty  )  is  written  as  f<pP^i  (2)  becomes 

(9)  PL^  /-  (/ -KpPn)*^ 

If  A/  is  changed  to  A/*  and  Xp  to  $  vrtiere  N*  m  cN  and  K*  ■ 
Kp/c  that  ■  NKp  »  the  equation  becomes 

(10) 

It  is  necessary  to  show  that  PL*  is  relatively  insensitive  to  changes 
in  the  factor  c  for  parameter  values  that  are  of  lntei«st. 
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It  p{x  t  y  )  is  small,  the  probability  of  no  hits,  {l  •  p  ,  may 
be  accurately  approximated  by  the  Poisson  form,  ,  in  which  the 

factor  c  cancels  out  completely.  Then  PL  becomes  approximately 
Even  if  P  is  not  small  enough  for  PL  to  be  independent  of  c  ,  the 
relative  change  in  PL  with  c  is  not  large  if  N  is  not  too  small.  Some 
sample  values  of  PL  for  various  values  of  c  are  as  follows: 


c  cN 

/  -  (/ -  p/c)^^( for 

Ns  /o) 

p  «  .025  .05 

.10 

.20 

.40 

1  10 

.224  .401 

.651 

.893 

.994 

2  20 

.222  .397 

.642 

.878 

.988 

4  4o 

.221  .395 

.637 

.872 

.985 

8  80 

.221  .39^ 

.63^ 

.868 

.983 

-pA/ 

/  -  e 

.221  .393 

.632 

.865 

.982 

Since  the  primary  objective  is  to  extend  the  Handbook  data  to  larger 
values  of.  N  t  it  appears  that  the  use  of  artificial  values  M*  and  Kp*' 
with  the  Handbook  graphs,  such  that  AJ^JC'p*  ■  AJKp  ,  is  a  reliable 
technique.  This  conclusion  is  supported,  as  described  earlier,  by  ex¬ 
amination  of  the  Handbook  graphs  and  the  results  of  special  machine 
computations. 

The  first  of  the  two  approximations  that  will  be  introduced  to  change 
the  form  of  the  basic  equation  is  the  substitution  of  the  Poisson  form, 

,  suggested  above  for  iX  •  p  )  in  (7)»  This  approximation  is  as¬ 
sumed  to  be  acceptable  if  N  is  larger  than  about  10  and  p  is  small. 

Here,  of  course,  p(A’,y)isa  variable  defined  by  (3),  but  it  is  al¬ 
most  always  small.  It  can  be  large  only  if  Kp  is  near  1  and  the  target 
dimensions  are  large  with  respect  to  the  ammunition  standard  deviations, 
as  shown  by  (8).  Since  the  second  approximation  will  require  that  the 
ammunition  standard  deviations  be  large  conqpared  to  the  target  dimensions, 
it  may  be  assumed  for  this  discussion  that  x,  y  )  will  be  small  even 
if  Kp  is  near  1. 

For  convenience  in  the  following  discussion,  (8)  is  written  as 

^  pi  X,  y)  =  ^p  (x)  It  (y)/2  TT 

I,ix)  "exp  d 4  (y)  =  fjexp  (-  f  v*)  dv 

Substitution  of  pix,  y)  from  (U)  into  (?)  «uid  replacement  of  (l 
by  e  then  give 

(12  )  PL»  i-  mpp  (-A/Kp  I,  4  /Z  Tt) 

Also,  before  introducing  the  second  approximation,  a  simple  change  of 
variables  is  convenient; 
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(13) 

In  terms  of  X  and  Y ^  the  original  expression,  (6),  becomes 

f  PL  exp (x^+  r^)]  dXcLY 

where  PL  is,  of  coarse,  expressed  in  terms  of  X  and  Y  ,  Further,  as 
suggested  by  the  comparison  of  curves  from  the  Handbook,  the  ratios  of 
ammunition  standard  deviations  to  delivery  standard  deviations  are  in¬ 
troduced  as  parameters  in  place  of  the  ammunition  stcuidard  deviations 
themselves ; 


(15) 

Then  the  limits  of  integration  in  (8)  may  be  written  in  terns  of  m  and 
n  and  the  new  variables  X  and  Y  as  follows : 

A  =  -A//W  -  XL/e^o^  C5f,  ^  ^  •  -y/'^  -  yt/zn 

a>=^x/rr»^rL/Zm<^  ,  '‘Y/n  Y^/zn 

For  the  second  approximation,  it  is  assumed  that  the  target  dimensions, 
XL  and  YLf  are  small  with  respect  to  the  ammunition  standard  deviations, 
and  •  The  integrals  in  (8)  may  then  be  removed  by  the  use  of 

I^i(u)cCu  ^  (0-A)f(i.A  -^-£3) 

Then,  in  terms  of  X  ,  Yi  m$  and  n,  the  Integrals  and  in  (ll) 
are  written 

s  (xl/^  cp;i^)exp 

it  =(yl/»  <^,)  e^p  (-  i  Yy„^) 

Substitution  of  these  into  (12)  gives 

(16 )  pl(x,  r)  ^  /-exp[-(^X^  r/zTTmn  c?? 

where  T  is  the  target  area  XL  YL  •  As  suggested  by  the  observations 
of  the  Handbook  data,  PL  can  now  be  written  in  terms  of  the  parameter 
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Z ,  where 


Z  «  /.  .5 

Also,  it  is  both  convenient  smd  reasonable,  as  discussed  previously,  to 
assume  that  /n*n.  The  common  ratio  can  then  be  denoted,  as  before,  by 
K4  ;  that  is. 


In  terms  of  and  A^,  then,  the  expression  for  PL  is 


(17)  PL  (x^  y)  =  /-  exp^-  C^/sTTA^^) exp  (-^ 

Since  this  is  the  only  place  where  any  parameters  are  involved  when  this 
expression  for  PL  is  substituted  into  (l^),  the  salvo  kill  probability, 

P  ,  may  be  regarded  as  a  function  of  just  Z  and  M  when  the  two  approxi¬ 
mations  are  valid.  The  first  approximation  is  valid  if  N  is  large 
enough  and  the  second  approximation  is  valid  if  XL/o^  and  YL/oi 
are  small  enou^.  These  are  essentially  the  same  conditions  that  wef^ 
derived  by  inspection  of  the  graphs  in  the  Handbook. 

When  the  approximate  expression  for  PL  from  (17)  is  substituted  into 
(li*’),  it  is  convenient  to  change  variables  to  polar  coordinates  (^,  5  ), 
where  ^  X  ^  +  Y^,  Integration  of  the  angular  variable  ^  then 
cancels  the  factor  l/27r  in  (l4),  giving 

(18 )  p  =  exp[j-  (z^37rM^)  exp  (-  py'z  exp  oL/i 

It  is  convenient  to  make  another  change  of  variable  and  two  more  changes 
in  notation; 

expC--^^^) ,  >s^z/s7r 

Then  the  approximate  expression  for  the  salvo  kill  probability  is 

/ 

(19)  /^ss  &xp  (-S  A  IV  ^JcCtV 

If  desired,  this  compact  form  of  the  basic  equation  could  be  used  in  a 
computation  program  to  find  values  of  P  for  cases  in  which  the  approxi¬ 
mations  are  valid.  There  is  actually  no  reason  for  doing  so,  since 
exact  values  of  P  have  been  computed  for  the  Handbook.  In  a  few 
special  cases,  however,  the  integral  in  (19)  assumes  familiar  forms  euid 
values  of  P  can  be  easily  calculated  without  machine  computation.  If 
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a  1,  for  example,  the  integral  becomes  a  simple  exponential  form, 
giving 

(20)  p(h=t)  ^  (j/s)(/-  e~^) 

If  >7=2,  the  integral  can  be  expressed  in  terms  of  the  standard  normal 
probability  function.  In  this  case. 


_/  . 

(21)  p(h=z)=  J  -  (-Tr/zs)^ 

-here 

Curves  of  P  as  a  function  of  Z.  corresponding  to  these  two  equations  are 
presented  in  Fig,  3;  values  for  a/)  a  1  were  computed  fran  (20 )  and  those 
for  M  =  0.707  from  (21 ).  Values  from  these  curves  agree  well  with 
those  from  the  Handbook  in  cases  where  the  approximations  might  be  ex¬ 
pected  to  be  applicable.  Tl'ie  curves  may  be  useful  to  the  reader  who 
wishes  to  interpolate  among  the  curves  of  Fig,  2. 

If  the  value  of  h  is  not  1  or  2,  the  integral  in  (19)  can  still  be 
evaluated  approximately  by  use  of  Sii^son's  rule,  although  it  may  be 
desirable  to  use  different  interval  sizes  in  different  parts  of  the 
range  of  integration.  If  values  of  P  for  Just  a  few  sets  of  values  of 
A/  and  Z.  are  desired,  the  amount  of  hand  computation  required  is  not 
prohibitive.  Values  computed  in  this  way  directly  from  (I9)  have  been 
used  to  plot  curve  a  of  Pig,  4,  which  shows  P  as  a  function  of  M  for 

Z  »  30. 

The  other  curves  in  Fig.  4  are  taken  from  the  H^dbook  and  are  all 
for  combinations  of  parameter  values  for  which  Z  ■  30.  The  curves  show 
the  effect  of  larger  targets  or  srasUJ.er  standard  deviations.  Curve  b, 
for  which  the  standard  deviations  are  rather  large  with  respect  to  the 
target  dimensions,  is  quite  close  to  curve  a,  the  graph  of  (19)  with 
Z  s  30,  but  the  other  curves  diverge  increasingly  from  a  and  b  at  the 
lower  values  of  M,  Uiis  is  to  be  expected,  since  if  the  deli  wry 
stsuidard  deviations  are  not  large  enough  with  respect  to  the  target 
dimensions  a  small  value  of  ^  meems  that  the  ammunition  steuidard  devia¬ 
tions  are  not  large  enough  for  the  second  approximation  to  be  useful. 

It  appears  from  Fig.  4  (and  adso  from  the  Handbook  graphs)  that  the 
effects  of  larger  targets  or  smaller  (telivery  steuidard  deviations  are  to 
increase  the  kill  probability  at  lower  values  of  the  ratio  Afcuid  to 
decrease  the  optimum  value  of  M  somewhat.  There  is  relatively  little 
change  in  the  mEucimum  value  of  kill  probability,  which  suggests  that  a 
graph  such  as  Fig.  1  should  be  reliable  except  when  the  delivery  standard 
deviations  are  quite  small  with  respect  to  the  target  dimensions.  Com- 
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parison  of  values  from  Fig.  1  auid  from  the  Handbook  tends  to  support 
this  conclusion. 

While  Fig.  1  was  plotted  from  Handbook  data,  presumably  it  is  also  a 
graph  of  the  function  of  Z.  obtained  by  maximizing  (19)  with  respect  to 
h  .  An  expression  for  this  rather  complicated  function  of  Z  has  not 
been  derived,  but  one  that  seems  to  approximate  it  rather  well  is 

(22)  =  [/- eyp  [- 


Actually,  this  expression  for  the  maximum  salvo  kill  probability  is  sug¬ 
gested  by  considering  a  somewhat  different  problem,  in  which  a  uniform 
distribution  of  borablets  over  an  area  rather  than  a  two-dimensional 
normal  distribution  is  assumed,-^  As  indicated  by  the  apparent  similarity 
of  the  functions  described  by  Fig.  1  and  by  (22),  the  maximum  salvo  kill 
probability  obtained  when  the  distribution  of  bomblets  is  optimized  often 
seems  to  be  quite  independent  of  the  particular  type  of  distribution 
assumed  for  the  borablets. 


^  Ihis  problem  is  described  in  Appendix  B  of  the  following:  U,  S, 
Naval  Ordnance  Test  Station,  Tactical  Use  and  Characteristics  of 
Effective  Free-Fall  Ordnanca  (u),  by  Roy  Dale  Cola,  M.E.  O'HelU,  M.R. 
O'Neill,  and  Richard  B,  Seeley.  China  Lake,  Calif.,  NOTS,  17  October 
1961.  (NAVWEPS  Report  78OI,  NOTS  TP  2797),  SECRET/NOFORN. 
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SUMMARY 

The  data  presented  in  NAVORD  REPORT  7019  (the  "Handbook" )  on  the  kill 
probability  of  cluster  weapons  against  unitary  targets  appear  to  include 
relationships  that  allow  the  ranges  of  values  of  parameters  to  be  ex¬ 
tended  beyond  the  Handbook  limits.  It  is  especially  in  the  data  for 
larger  numbers  of  bomblets  and  larger  ratios  of  standard  deviations  to 
target  dimensions  that  the  relationships  are  observed. 

The  salvo  kill  probability,  P  ,  may  be  regarded  as  a  function  of  the 
product  NKp  ,  if  A/  is  larger  than  about  10,  (  A/  is  the  number  of  bomb- 
lets  and  Kp  is  the  probability  that  a  bomblet  hit  is  a  kill. )  TSiis  means 
that  if  hi  is  beyond  the  range  of  Handbook  values,  P  may  be  determined 
by  using  the  Handbook  graphs  with  "artificial"  values  hi*  and  /a*8uch 
that  hJ^’K'p^.:.  NKp  . 

Further,  if  the  ammunition  standard  deviations,  suid  <7^  ,  are 

each  larger  than  about  half  the  corresponding  target  dimension^  the 
salvo  kill  probability  may  be  regarded  as  a  function  of  Just  two  para¬ 
meters  Z.  and  M  ,  where 

Z  s  /.S  NKp  t/(7p^  <JI^/C7p^  - 

(Hie  factor  1.5  in  Z  corresponds  to  the  fixed  ratio  of  Op^  to 
used  in  the  Handbook,  and  ^  are  the  delivery  standard  devia¬ 

tions  and  T  is  the  target  area.)  Interpolation  or  extrapolation  of 
the  Handbook  data  is  thus  easily  accomplished  by  choosing  artificial 
values  that  leave  Z.  and  M  unchanged  but  that  are  within  the  Hand¬ 
book's  ranges  of  variables. 

This  approach  works  especially  well  if  values  of  the  maximun  se^vo 
kill  probability,  P,  ,  are  desired.  Also,  since  ^  is  a  function  of  z 
only,  the  functional  relationship  can  be  shown  on  a  graph,  such  as  Flg.l 
of  this  report.  Values  of  P^  from  this  graph  appear  to  be  accurate 
within  about  0.01  if  N  is  greater  than  about  10  and  and  axe 

greater  than  about  half  the  corresponding  target  dimensions. 

A  graph,  such  as  Fig.  2,  can  also  be  constructed  to  show  P  as  a 
function  of  the  two  variables,  Z  and  A/,  Values  of  P  from  this 
graph  seem  to  be  accurate  within  a  few  hundredths  if  the  ammunition 
dispersions  are  such  that  either  A/  is  not  less  than  optimun  or 
and  are  greater  than  about  half  the  corresponding  target  dimensions. 
The  io^rtance  of  using  optlmin  values  of  M  (that  is,  of  incorporating 
optimiia  ammunition  dispersions  into  the  weapon)  is  also  shown  by  these 
curves.  If  a  cluster  weapon  is  to  have  a  relatively  hic^  salvo  kill 
probability  the  desired  asmunltlon  standard  deviations  should  be  about 
as  large  as  the  expected  delivery  standard  deviations. 
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These  empirically  derived  techniques  are  generally  supported  when 
the  basic  equation  used  in  the  preparation  of  the  Handbook  Is  examined. 
Two  approximations  that  might  be  expected  to  hold  if  the  number  of  bomb- 
lets  is  large  enough  and  the  ratios  of  target  dimensions  to  standard 
deviations  are  small  enouj^  result  in  a  rather  compact  form  of  the 
basic  equation. 

It  is  hoped  that  the  methods  described  in  this  report  will  be  help¬ 
ful  when  the  effectiveness  of  salvo-delivered  and  clustered-bomblet 
weapons  is  to  be  estimated.  It  is  possible  to  take  account  of  rather 
large  variations  in  parameter  values  without  requiring  additional  long 
(costly)  use  of  high-speed  computing  equipment. 
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FIG.  1.  Maximum  Kill  Probability  Vs,  Z,  (Optimum  Ammunition  Dispersion) 


•  KILL  PROBABILITY 


1.0 


■«|g8— ” 

i^^raBF 


[2  =  20 


2  =  15 


Volues  of  Sr**  Probability  Vs.  the  Ratio  m 


17 


NAVWEPS  REPORT  8075 


18 


FIG*  3*  Two  Curves  of  Salvo  Kill  Probability  Vs,  Z  Computed  Directly  from  the  Approximate  Equati 
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Code  466  (1) 

1  Fleet  Anti-Air  Warfare  Training  Center,  San  Diego 
1  Naval  Air  Development  Center,  Johnsville 
1  Naval  Missile  Center,  Point  Nugu  (Technical  Library) 

1  Naval  Ordnance  Laboratory,  Corona 

1  Naval  Ordnance  Laboratory,  White  Oak  (Technical  Evaluation  Group) 

1  Naval  Postgraduate  School,  Monterey 
1  Naval  Research  Laboratory 

1  Naval  Weapons  Evaluation  Facility,  Kirtland  Air  Force  Base 

1  Naval  Weapons  Laboratory,  Dahlgren  (Technical  Library) 

2  Naval  Weapons  Services  Office 

1  Operational  Test  and  Evaluation  Force 

1  Bureau  of  Naval  Weapons  Fleet  Readiness  Representative  Pacific,  Naval 
Air  Station,  North  Island 

3  Chief  of  Ordnance 
ORDFA  (1) 

ORDTB  (1) 

ORDTS  (1) 

2  Aberdeen  Proving  Ground  (Ballistic  Research  Laboratories) 

1  Weapon  Systems  Laboratory 

6  Army  Missile  Command,  Redstone  Arsenal 
ORDDW-IDE  (1) 

Rocket  Development  Division,  Test  and  Evaluation  Branch  (1) 

Redstone  Scientific  Information  Center  (4) 
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1  Diamond  Ordnance  Fuze  Laboratories  (Ordnance  Development  Laboratory) 

1  Picatinny  Arsenal  (Library) 

1  White  Sands  Missile  Range  (ORDBS-Technical  Library) 

2  Headquarters,  U.  S.  Air  Force 
AFDRD-AN  (1) 

AFDRD-CC  (1) 

1  Strategic  Air  Command,  Offutt  Air  Force  Base 

2  Air  Proving  Ground  Center,  Eglin  Air  Force  Base 
P6APZ,  Technical  Library  (1) 

WWTA,  Weapons  and  Target  Director,  Detachment  4  (1) 

1  Holloman  Air  Force  Base 

1  Space  Systems  Division,  Air  Force  Systems  Command,  Los  Angeles 
1  Assistant  Secretary  of  Defense  (R  &  E)  (Technical  Advisory  Panel  on 
Ordnance ) 

1  Defense  Atomic  Support  AgSncy,  Ssndia  Base 
1  Weapons  Systems  Evaluation  Group 
10  Armed  Services  Technical  Information  Agency  (TIPCR) 

1  Langley  Research  Center  (Library) 

1  Aerojet-General  Corporation,  Azusa,  Calif.,  via  BuWepsRep 
1  Applied  Physics  Laboratory,  JHU,  Silver  Spring 

1  Armour  Research  Foundation,  Chicago  (Document  Librarian  for  Department  M) 

1  Autonetics,  A  Division  of  North  American  Aviation,  Inc.,  Downey, 

Calif.  (Department  93,  F.  H.  Gardner) 

2  Convair,  San  Diego  (Engineering  Library) 

1  Douglas  Aircraft  Company,  Inc.,  Long  Beach  (Technical  Library) 

1  Douglas  Aircraft  Company,  Inc.,  Santa  Monica,  Calif.  (Missile  Division) 

1  Jet  Propulsion  Laboratory,  CIT,  Pasadena  (Dr.  W.  H.  Pickering) 

1  Lockheed  Aircraft  Corporation,  Burbank,  Calif. 

1  Lockheed  Aircraft  Corporation,  Missile  Systems  Division,  Van  Nuys,  Calif. 
1  Los  Alamos  Scientific  Laboratory  (Reports  Librarian) 

1  New  Mexico  Institute  of  Mining  and  Technology,  Socorro  (Dr.  Workman) 

1  North  American  Aviation,  Inc.,  Los  Angeles  (Weapons  Systems  Group) 

1  North  American  Aviation,  Inc.,  Space  (  Information  Systems  Division, 
Downey,  Calif.  (Department  56,  Group  55) 

1  Operations  Research  Office,  JHU,  Bethesda,  Md. 

1  Project  THOR,  JHU,  Silver  Spring 

1  Stanford  Research  Institute,  Poulter  Laboratories,  Menlo  Park,  Calif. 

1  The  Boeing  Company,  Seattle  (Structural  Staff  Engineer) 

1  The  Rand  Corporation,  Santa  Monica,  Calif.  (Aero-Astronautics  Department) 
1  University  of  Chicago,  Institute  for  Air  Weapons  Research,  Chicago 
1  University  of  Denver,  Denver  Research  Institute,  Denver  (A.  M.  Krill) 
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